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Abstract. We study the properties of a heterogeneous, chiral granular rotor that is 
capable of performing useful work when immersed in a bath of thermalized particles. 
The dynamics can be obtained in general from a numerical solution of the Boltzmann- 
Lorentz equation. We show that a mechanical approach gives the exact mean angular 
velocity in the limit of an infinitely massive rotor. We examine the dependence 
of the mean angular velocity on the coefficients of restitution of the two materials 
composing the motor. We compute the power and efficiency and compare with 
numerical simulations. We also perform a realistic numerical simulation of a granular 
rotor which shows that the presence! of non uniformity of the bath density within 
the region where the motor rotates, and that the ratchet effect is slightly weakened, 
but qualitatively sustained. Finally we discuss the results in connection with recent 
experiments. 



PACS numbers: 05.20.-y,51.10+y,44.90+c 
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1. Introduction 

Much is known about the properties of Brownian ratchets and their apphcations in 
physics and biology. Smoluchowski's [Ij original proposal consisted of a device with 
four vanes connected to a ratchet and pawl. When immersed in a thermalized bath of 
particles the device can apparently rectify the thermal fluctuations to produce useful 
work. Feynmann[2], however, demonstrated that this is not possible if the entire system 
is at thermal equilibrium. The device only performs work if the vanes and the ratchet and 
pawl are maintained at different temperatures, Ti and T2, respectively with Ti > T2. 
It is now understood that the general requirements for a functioning motor are the 
absence of both time and spatial symmetry. The former occurs if there is a breakdown 
of detailed balance and the latter may be due to the intrinsic asymmetry of the object 
itself or to an external force [3], HJ [5] . 

In the last few years, several proposals for granular motors have appeared. The 
dissipative nature of inelastic collisions leads to an automatic breakdown of time reversal 
symmetry and the spatial asymmetry can be created in various ways. In 2007 Cleuren 
and Van den Broeck [H] and independently, Costantini, Marconi and Puglisi [71 [S] 
proposed the same model granular motor, viz. an isosceles triangle composed of a 
homogeneous inelastic material constrained to move along a line. When immersed in a 
bath of thermalized particles, the undirected fluctuations of the bath particles induce a 
directed motion of the triangle. The drift velocity is proportional to (1 — a) ^m/M 
where a is the coefficient of restitution characterizing the collision between the tagged 
particle (of mass M) and the bath particles (of mass m). This result shows that the 
phenomenon is specific to granular particles, because the drift velocity vanishes when 
a = 1. Also, because the granular temperatures of the tagged particle and the bath 
particles are comparable, whatever the mass ratio, < {V— < V >Y > decreases as 
m/M (for a given bath temperature), and this results in a signal-to-noise ratio that 
vanishes for large values of ^Jm|M. 

Cleuren and Eichhorn [9] later proposed an alternative model consisting of a 
homogeneous rotor that is constrained to rotate about an off-center axis. If the 
material composing the rotor is inelastic, a net rotation is obtained and it displays the 
same dependence on a and m/M a,s the triangle, which complicates the experimental 
observation of this phenomena. 

Fortunately, this difficulty can be overcome by using a heterogeneous device 
constructed from two different materials with different coefficients of restitution. This 
leads to a strong ratchet effect [lU]. Costantini et al. [TT] proposed possibly the 
simplest model of a granular motor of this type, the asymmetric piston composed of 
two materials with different inelasticities. Starting from the Boltzmann equation, they 
proposed a phenomenological approach based on the evolution of the first three moments 
of the velocity distribution. Although their theory is in good agreement with numerical 
simulation results in some cases, it breaks down in the limit of large piston mass. 

Recently, Eshuis et al. [12] succeeded in constructing a macroscopic rotational 
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ratchet consisting of four vanes that rotates in a granular gas. They performed two 
kinds of experiments: in the first, the vanes were symmetric and no net angular velocity 
was observed. In the second, asymmetry was induced by coating each side of the 
vanes differently. For sufficiently large granular temperatures of the bath particles, 
a net rotation is observed. We will discuss the implications of this experiment in the 
concluding section. 

In a recent article [T3j we presented a mechanical approach that gives the average 
force acting on a heterogeneous particle moving at a given velocity. The mechanical 
approach is consistent with the Boltzmann equation and in the Brownian limit, m/M — >■ 
0, it leads directly to an exact expression for steady state drift velocity. 

From both technological and fundamental perspectives, the power and efficiency of 
thermal engines are of major interest. Indeed, Carnot's analysis of an idealized heat 
engine is a landmark event in the history of thermodynamics. His famous result for the 
efficiency depends only the temperatures of the heat reservoirs, but applies strictly to a 
reversible - and hence quasi-static - process. The efficiency at non-zero power was first 
considered by Curzon and Ahlborn[l3] and later generalized by van den Broeck [15]. 
The efficiency of Brownian motors has also been discussed [T7] , but there has been 
little discussion of the power and efficiency of granular motors [13] . Such considerations 
are, nevertheless, of great importance for the practical realization of these devices. 

It is the purpose of this article to apply the approach first proposed in [13] to 
the chiral rotor. We also investigate the effect of bath inhomogeneities with a realistic 
numerical experiment. Finally, we show that the mechanical approach provides an 
explanation for the presence of a threshold for the rotation of the vanes in the experiment 
of Eshuis et al. [12]. 

2. Model 

The chiral rotor is composed of two materials with coefficients of restitution a+ and 
a_: See figure [Tj Collisions of the bath particles with the former result in a positive 
torque (anticlockwise sense), while collisions with the material with a_ result in a 
negative torque. The device is immersed in a two-dimensional granular gas composed 
of structureless particles at a density p each of mass m. 

An important assumption of the model is that collisions between the motor and 
the bath particles do not modify the velocity distribution of the latter. This neglect of 
recoUisions is consistent with the foundations of Boltzmann's kinetic theory. We let v± 
and f II denote the components of the gas particle's velocity perpendicular and parallel to 
the surface of the motor, respectively. The granular temperature of the bath is defined 
as Tb = m < f 5^ >= m < v"^^ >. Since v\\ is irrelevant in the Boltzmann-Lorentz 
description of the rotor [18j, for sake of simplicity we denote v± by v in the rest of the 
section. 

Wherever possible, we will give results for an arbitrary bath particle velocity 
distribution, 0(t>) = J (it>|j0(v). Sometimes, though, it will be necessary to assume 
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Figure 1. The chiral rotor is constructed from two different materials with coefficients 
of restitution a+ (red) and a- (blue). When immersed in a bath of thermalized 
particles, it acquires a net rotation around its axis. 



a particular form and in this case we will use a Maxwell-Boltzmann or Gaussian 
distribution: 



1 / V 



,2 



(f)MB{v) = exp -— ^ , (1) 



where Vth = \/TB/ra. While granular gases often have non-Maxwellian velocity 
distributions, it is possible to devise experimental situations where the Gaussian 
distribution is observed [121 120] • 
The collision equations are 

"'^ -^(..-n.)f I (2) 

1 + mx^ \ —I I 

where Q is the angular velocity, / is the moment of inertia and —L/2 < x < L/2 
is the algebraic distance of impact from the center (for simplicity, we neglect collisions 
between bath particles and the caps of the rotor). The granular temperature of the rotor 
is given by Tg = / < {Q— < Q >)^ >. Because the mass of any section of thickness 
Ax perpendicular to the long axis of the rotor is constant, the moment of inertia is the 
same as for a homogeneous rotor, / = ML^/12, where M is the total mass. The energy 
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loss on collision is 

ml{l-a^){v-nxf 

= 2(/ + mx^) 

The reconstituting velocities are 

V I I + mx^ \ —I I 

3. Boltzmann-Lorentz Equation 

The time dependent Boltzmann-Lorentz equation for the chiral rotor is 
a /.L/2 ^-oo r f(Q**-f) f(Q* 

— /((];t) =p dx dv\v-xn\ eiv-xn y^ ' ' (t){v**) + eixvt - vY-^. 

/L/2 POO 
dx / dv \v — x^l\(j){v) 
■L/2 J-oo 

where 6{x) is the Heaviside function. By introducing the variable y = {v — Qx)/a^, 
y = {Qx — v)/a- and exploiting the symmetry of the object we can write the equation 
in the more explicit and useful form: 



2pdt 



f{il;t) = dx dyyfiil-— -mxy ] (p [ ilx + —— 

Jo Jq \ I + mx^ J \ I + mx^ 

+ / dx dyyf{il + — -mxy ]<p\ilx — —y 



JO 

L/2 POO 



fl^l^ fOO 

-fin) dx dyy{(l){nx + y) + (t>{nx-y)) (6) 
Jo Jo 

When the bath distribution is Gaussian, equation [TJ one can confirm that for 
a rotor with a+ = a_ = 1, the solution in the steady state is given by /(i^) = 
^I/2TiTBeM-m''/2TB). 

3.1. Numerical solution 

The Direct Simulation Monte Carlo (DSMC) method [21] is often used to obtain 
numerical solutions of the Boltzmann equation and has been applied to granular systems 
[22] . An alternative is the Gillespie method [23]. In this application, the linear nature 
of the Boltzmann-Lorentz equation makes it well suited to solution by iteration. To 
achieve this, we rewrite the equation in the following form: 

"^''^ , , .(^\ ( ^ 1 + a+ \ . ( ^ I - a+mx^ 



mx"^ 



r-L/2 POO / \ -\- a \ f 

+ dx dy yf^"^^ f H —mxy ) f fix 

Jo Jo V I + mx^ J V 



/ + mx^ 
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where Z*^"-* denotes the nth iteration and 

/•L/2 POO 

iy{n)= / dx dy y{(j){VLx + y) + (j){VLx - y)) (8) 
Jo Jo 

so that 2pi>{Q) is the colhsion rate when the motor rotates with angular velocity Q. We 
take as the initial guess, /*^°\ either a Gaussian distribution or the previous converged 
solution for different parameters. The error is defined as 

oo 

_/W(())]2^n (9) 

'OO 

In the results reported here we took the solution as converged when e < 10"^^. 

We show some results in figure [2j As the rotor becomes heavier, the distribution 
becomes more Gaussian with the maximum at the drift velocity Q*. Conversely a light 
rotor has a markedly non-Gaussian distribution that becomes more and more skewed 
with decreasing mass. For a sufficiently light rotor, the most probable angular velocity 
is opposite in sign to the mean value. This is similar to the behavior observed for the 
asymmetric piston [7j. 

The increasingly non-Gaussian behavior can be quantified by examining the 
skewness and kurtosis, k, of the distribution: See figure [3] (for convenience the excess 
kurtois, defined as k — 3, is shown). While the excess kurtosis is nearly zero beyond 
I/mL"^ = 1, there is some residual skewness for large rotor masses. We also note the 
the granular temperature of the rotor varies little beyond / / mL^ = 1 . 

3.2. Kramers-Moyal expansion 

Henceforth we focus on the Brownian limit as it is more interesting for the experimental 
application. As with the piston [TH] we can perform a series expansion of Boltzmann 

equation by introducing the small parameter e = the variable z = {Q — Q*) / e 

where Q* is the drift velocity in the Brownian limit, and a rescaled velocity distribution 
F{z). At first order we obtain 

POO /'L/2 

dy dxxy^{{l + a+)(f){xn* + y)-{l + a-)(l){xn* -y)) = 0, (10) 

Jo Jo 

which can be solved numerically to obtain Q*. At second order we obtain a Fokker- 
Planck equation for the rotor angular velocity distribution that gives, in the original 
variables 

/(^])oce-^(^-^*)'/2T,_ 
where Tg/m is given by 

7, ^ jl^'^dx j;^dyx'y\{l + a.)mn*x-y) + {^ + a^)mn*x+y)) 

^ 4:J^^^dx J^dyx^y{{l + a^)(l){n*x-y) + {l + a+)(l){n*x + y)) 

For a rotor with a+ = = in a Gaussian bath, equation ([T]), we find 

Tg/Ts = 0.72960 which is close to the value for the corresponding asymmetric piston 
{Tg/TB = 0.73460). 



Kinetic analysis of a chiral granular motor 



7 



0.5 r 




0.2 r 



iG 0.1 




Figure 2. Angular velocity distributions for a chiral rotor with q:_ = 0, a+ = 1 
obtained by iteration of the Boltzmann equation. Top: 1 < / < 10 in steps of 1; 
Bottom 0.1 </ < 1.0 in steps of 0.1 



4. Torque-based approach 

The instantaneous impulse exerted on the rotor of angular velocity Qhj a. collision with 
a bath particle of velocity w on a face with is 



/_ = m' - m 

Imx 



I + mx"^ 

Similarly for a collision on a face with a 

Imx 



(/ + a^){Qx — v), V < Qx 



{I + a-^-){v — ^Ix), ^Ix < V 



(13) 



(14) 



Assuming that successive collisions between bath particles and the granular motor are 
uncorrelated, averaging on different collisions at a given Q consists of integrating the 
impulse times the rate that the specified face collides with a particle moving with a 
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Figure 3. First moment (second from top, right hand side), second moment (third 
from bottom, right hand side), skewness (second from bottom, right hand side) and 
excess kurtosis (bottom curve) of the velocity distribution for = 0, Q!-|_ = 1 as a 
function of /i. The dashed Hne shows the granular temperature of the rotor. 



velocity v. Evaluation of the torque requires an integral over x, in addition to the bath 
particle velocity distribution, i.e. 

L/2 i>nx 

dx I dv I_{n,v){nx-v)(j){v) (15) 



T_{n) = 2p 



L/2 



dx 



■oo 
oo 



dv /+(f2, v){v — Qx)(f){v) 



(16) 



Ha; 



We suppose that, in addition to the fluctuating force resulting from collisions with 
the bath particles, the rotor is subject to a constant external torque. Text- The average 
of the net torque performed over all possible angular velocities of the rotor is equal to 
zero in the stationary state. Since the successive colhsions are assumed uncorrelated, 
this condition can be expressed as, 







(17) 



-Xn)>= / dnf{n)T{n) + T,, 

J — oo 

where = T^{Q) + r+(r2). It is convenient to rewrite this as 

<r„et(^^) >=^(<^] >) + !/ + Text = (is) 

where r(< Q >) is the mean torque when the rotor has constant angular velocity, < Q > 
and Tf is the mean torque resulting from the fluctuations of the angular velocity around 
the mean value, < Q >. Now by expanding T{Q) about Q*, where r(f2*) = we have 
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Figure 4. Dimensionless mean angular velocity as a function of the mass ratio 
M/m for a+ = 1 and «_ = 0. The joined circles show the simulation results, the 
solid horizontal line is the exact result in the Brownian limit (the solution of r(ri) ~ 



with r(r2) given by equation (26)) and the dashed line shows the first order correction, 



equation (23) 



that 



+ -(Q - Q*fT"'(Q*) + ... > -r(< Q >) 
6 



So to lowest order the angular velocity is given by 

<n>=n 



^ ir"(fi*)<(fi-fi*)2>+rext 



Since I < {Q — Q*y >^ Tg, we have that in the Brownian limit, / — )■ cxo, 

<fi>=^]*-^ext/^'(fi*) + o(l//) 



(19) 



(20) 



(21) 



Making the substitutions y = Qx — v and y = v — Qx, we find the following 
expression for the total torque in the Brownian limit, mL"^ / 1 — t- 0: 

i'L/2 POO 

T{VL) = 2mpi dx dyxy^{{l + a+)(l){xn + y)-{l + a.)(f){xn-y)) (22) 
Jo Jo 

We obtain the mean angular velocity, Q*, by setting the torque equal to zero: r(i7*) = 0. 

The result, as well as simulations for finite ratios are shown in figure |4j Note that, 
unlike the piston [13], we do not observe a maximum as the moment of inertia decreases. 

We can also calculate corrections to the Brownian limit by using equation ( [20) ) 

1 T"{n*) Tg 



2 r{Q*) I 



(23) 
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Figure 5. Dimensionless torque as a function of the mean angular velocity fl for 
Q!_|_ — 1 and a_ = 0. The torque is positive for 51 = and vanishes for a positive value 
of il, which corresponds to the solution of r(ri) = with T{il) given by equation (26 1. 
At large angular velocity, the torque is opposite in sign to the angular velocity. 



Tg, T'{Q*) and r"(17*) are functions of a+ and a;_. Figure |4] shows that, for a;_ = 0, = 
1, the first order correction provides a reasonably accurate description for M/m > 5, 
but for smaller values it diverges in contrast to the exact result that approaches a finite 
value. 

For an arbitrary symmetric bath distribution, the torque on a stationary rotor in 
the Brownian limit is 

T{0) = lmpL'^{a+-a.) <v^ > (24) 
8 

In an experimental situation, this must exceed the solid friction for the motor to function. 

For small variations of the angular velocity, the torque exerced by the particle bath 
is proportional to the deviations from the steady state T{Q) = {Q* — Q)r'{Q*) where 
r'(f2*) is always a nonzero negative constant. 

At large angular velocities, one obtains that 

m) ^ ^ ^(1 + ^t)p^^ ^2^ ^ ^ ±^ (25) 
32 

The damping is quadratic at large values of the velocity (the piston has a similar 
behavior) . 

For a Gaussian distribution in the Brownian limit, we can obtain an analytical 
expression for the torque: 



—FYT- = + 8) - (2 + a+ + a_) 
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l(f]2 ^ 8 - i^)er/(-^) + - A ™ + i) exp(-fiV8) 



(26) 



where f2 = VLL/vth (see figure |5j). Expanding in powers of Vt to obtain 

r 1 .1/2" 



from which we can obtain an analytical estimate for the angular velocity in the steady 
state 

f2* = ^./^ (28) 
2V22 + a+ + a_ ^ ^ 

We note that this differs from the result for the asymmetric piston by a constant 
factor. This result also shows that, for given coefficients of restitution and bath 
properties, the rotational velocity is inversely proportional to the length of the rotor. 

4-1 ■ Dependence of the mean angular velocity on the coefficients of restitution 

We examine the mean angular velocity in the Brownian limit as a function of the 
coefficients of restitution. It is easy to show starting from r(r2*) = with T{VL) given by 



equation (22) that fi*(Q;+ = x,a^ = y) = — r2*(a;+ = y,a_ = x), regardless of the bath 
distribution. Of course this result is expected from simple symmetry considerations. In 
figure |6] we show fi*(a+ = x, a_ = 0) and f2*(a+ = !,«_.= x) as a function of x in 



the case of a Gaussian bath distribution. The former increases as x increases, while 



the latter decreases. They are equal for a particular value of x and from equation (26) 
we may show that this occurs for x = \/2 — 1. These results illustrate that the mean 
angular velocity is not simply proportional to the difference in coefficients of restitution. 
Specifically, we can construct two rotors. The first is composed of an elastic material 
(a;+ = 1), and another material with a coefficient of restitution of x = \/2 — 1 or 
a difference of 2 — \/2 = 0.5858. This has the same drift velocity as a second rotor 
composed of a completely inelastic material, (a;_ = 0), and one with a_ = x or a 
difference of 0.4142. 

These functions are very nearly symmetric i.e., to a good approximation they satisfy 
n*{a+ = X, a_ = 0) + n*{a+ = 1, a_ = x) = fi*(l, 0) (29) 
Figure |6] also shows that the exact solutions are well-approximated by the second 



order estimate equation (28) 



4-2. Power and efficiency 

A defining characteristic of granular motors is their ability to extract energy from 
the bath in the form of mechanical work. We now proceed to the calculation of the 
mechanical work, as well as the dissipated power due to the inelastic collisions between 
the granular motor and the bath particles. The ratio of these two quantities defines the 
efficiency of the motor. 
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Figure 6. il*{a+ — x,a^ — 0) and il*(Q:+ — 1,q:_ — x) and their sum as a function 



of X. The sohd lines show the solution of r(f2*) = with T(D,) given by equation (26) 



and the dashed lines show the approximate solution, equation (28) 



If the external torque in equation (17) is non-zero, work is being done by, or on, 
the rotor. The power is given by 

W = Text < > (30) 



In the Brownian limit we can estimate the average drift velocity using equation (21). 
Alternatively, we can take < Q > as the independent variable: 

W = T{<n>)<n> (31) 



with r(f^) given by equation (22 ). Of course, this approach is only valid in the Brownian 
limit. In general Vr(n) has a parabolic form and is maximum for < < When the 
angular velocity is greater than Q* the power is negative, implying that it is necessary 
to drive the rotor externally in order to maintain the motion. 
At second order 

)Q- ^{2 + a+ + a^)J-n^ (32) 



mpLv'^f^ 8^^ 12^ ' ^' ' y 71 

Figure [7] shows that this provides a good approximation. The maximum power occurs 
at Qrnp = ^* /2 and is given by 

mpLv'ff^ 128 2 + a+ + a 
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Figure 7. Power as a function of the rotor angular velocity for a system with a_ = 
in the Brownian limit. The solid and dashed lines show the exact (equation (26) in 
equation (31)) and approximate (equation (32)) values, respectively, a^^ = 1,0.7,0.5 
in order of decreasing amplitude 



From equation ([s]), we find that the rate of energy dissipation resulting from 
colhsions between the bath particles and the rotor is 

/•L/2 POO 

Etotin) = -mp dx dy y^'lil - a^)(l)inx - y) + (1 - al)(f)inx + y)] (34) 
Jo Jo 

The efficiency, or the fraction of the dissipated energy that is converted into work, 

is: 

W 

ri{a+,a.,n) = (35) 

For an arbitrary bath distribution in the Brownian limit we find that 
r]{n, a+, a_) _ J^^"^ dx dy xy'^{{l + a+)0(a;fi + ?/)-(! + a_)(t){xn - y)) 



2^ /,^/' dx dy y3((i _ al)<t>ixn + y) + il- aimxn - y)) 



(36) 



/o ""^ Jo 

For a Gaussian bath distribution this has a maximum value of 4.3335% for = 1, a_ = 
0,f2 = 0.29864. Figure [s] shows the efficiency as a function of the dimensionless mean 
angular velocity obtained by simulation for several mass ratios, M/m as well as the 



theoretical result in the Brownian limit, equation (36) 
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Figure 8. The efficiency as a function of the dimensionless mean drift angular 
velocity for a system with = 0,a+ — 1. The solid curves show simulation results 
for M/m = 1,2,5,10 (bottom to top), while the dashed line shows the theoretical 
prediction in the Brownian limit. The bath particle velocity distribution is given by 
equation Q. 



This mechanical analysis indicates that heterogeneous granular particles are able to 
produce significant noise rectification. Provided that the solid friction about the axis is 
not too large, the chiral rotor should be a good candidate for an experimental realization 
of a Brownian granular motor. 



5. Generalization to an n-vaned rotor 



If one retains the assumption that the bath remains homogeneous and is not influenced 
by the presence of the rotor, it is easy to generalize the above analysis to a rotor 
composed of n vanes. In the Boltzmann-Lorentz equation (|6]), the factor of l/2p on 
the left hand side is replaced by 1/np. In the steady state the entire left hand side is 
equal to zero, so the steady state angular velocity is unchanged. In the expression for 
the torque, the factor of 2 in equation (22) is replaced by n. But since this factor does 
not affect the solution of r(fi*) = 0, the number of vanes has no effect on Q*. We note, 
however, that the assumption of bath homogeneity certainly deteriorates as the number 
of vanes increases. 
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Figure 9. Top view of the simulation cell of size 20 x 20 sphere diameters: The 
spherocylinder has a radius R = 2/3r (r is the radius of the spheres) and a total 
length L + 2R where L = Hi?. 

6. A numerical experiment 

The analysis presented in the preceeding sections assumes that the bath is perfectly 
homogeneous. To investigate effect of relaxing this assumption, we now consider a 
realistic simulation of a quasi two-dimensional rotor. Burdeau et a/. [21] developed a 
model that accurately reproduces an experimental study of a vibrated bilayer. The 
first layer consists of densely packed heavy granular particles while the second layer is 
composed of light particles with an intermediate coverage. Both experimentally and 
in the simulation, one observes a horizontal velocity distribution function of the light 
particles that is very close to a Gaussian 

To this vibrating bilayer system, we have added here a three dimensional 
spherocylinder constrained to rotate around a vertical axis: See figure [9j The 
spherocylinder cannot move horizontally or vertically: the rotation around the Oz axis 
is the sole degree of freedom. Figure [TO] is a projection of the spherocylinder together 
with the geometrical parameters. 





L/2 






/\ Rs 


O 



















Figure 10. Geometry of the spherocylinder and a sphere during a collision 
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The simulation uses a Discrete Element Method where collisions between particles, 
collisions between spherocylinder and particles as well as collisions with the vibrating 
plate are all inelastic. In addition, the particles are subject to a constant acceleration 
due to the (vertical) gravitational field. The visco-elastic forces are modeled by the 
spring-dashpot model 

The systems consists of A^^i spheres placed in the first layer and N2 spheres in the 
second layer. By choosing the system parameters as close as possible to the original 
experimental system[T9], no mixing is observed between layers and the horizontal 
velocity distribution of the second layer is very close to a Gaussian. The force along 
the line joining the two centers that dissipation into account is a damped harmonic 
oscillator defined as: 

Fn = -Ki - Ini (37) 

where ^ is the overlap between two particles, A;„ is related to the stiffness of the material, 
and 7„ to the dissipation. This force model allows one to very easily tune certain 
quantities in the simulation, namely the normal coefficient of restitution e„ (which is 
then constant for all collisions at all velocities in our system). Relevant values of Cn 
and tn determine the values of and 7„, independently of the velocities. The collision 
duration, t„, provides a microscopic characteristic time. The simulation time step is 
taken as At = 10~^s and the mean duration of a collision as tn = 10~^s. A frictional 
force between spherical particles is added by taking the tangential component of the 
force as Ft = — min(|ktC|, |/^Fn|) where kt is related to the tangential elasticity and ( 
is the tangential displacement when the contact was first established. We used a ratio 
kt/kn = 2/7 with ^ = 0.25. In order that the spherocylinder only undergoes collisions 
with particles of the second layer, we choose the vertical height of its center of gravity 
a.s h = 3r + R. 

We first performed simulations for a homogeneous spherocylinder, where no mean 
angular velocity appears. If the linear dimension of of the simulation cell is at least 
twice the length of the spherocylinder, the boundary conditions do not infiuence the 
kinetic properties of the latter. We then considered a chiral spherocylinder made of two 
materials of coefficient of restitution a+ = 0.8 and a- = 0.2. 

The motor effect is clearly evidenced in this numerical experiment: See figure [TT] 
The Boltzmann-Lorentz approach can be extended to objects of finite area (see Ref. [25]). 
In the limit of large L compared to the bath particle radius and the radius of the cap. 



the mean angular velocity varies as 1/L as shown in figure 11 We also monitored 
the density of the bath particles in the region of the spherocylinder. Figure 12 shows 
the near exclusion of particle centers in the range < r < R + Rg due to the steric 
effect (the fact that the density is not strictly zero is due to out-of-plane collisions 
resulting from the three dimensional nature of the system). The vertical dotted line 
corresponds to r = R + Rg. One observes a small depletion in the region swept out by 
the rotating spherocylinder followed by an increase to a maximum at a distance close 
to {L + Rs)/2 + R = 9.5Rs i.e., a distance where bath particles collide with the ends 
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of the sphero cylinder. Beyond this distance the density rapidly approached the bulk 
value. These effects are similar for homogeneous and heterogeneous spherocylinders. For 
comparison the straight dashed line corresponds to the assumption of a homogeneous 
bath used in the Boltzmann-Lorentz description. 




Figure 11. Mean angular velocity versus the length of the spherocylindcr with a+ = 1, 
a_ = 0.25, and spherocylinder mass Ms = 40m. 



7. Conclusion 

We have applied kinetic theory and realistic numerical simulation to study the properties 
of a chiral heterogeneous granular motor. As with the previously studied piston, a 
strong ratchet effect may be produced under certain conditions. We also developed a 
mechanical analysis that provides an exact description in the Brownian limit where the 
motor is much more massive than a bath particle. 

The mechanical approach may also be used to assess the robustness of the ratchet 
effect. In an experimental situation, a dry friction force between the rotor axle and 
its bearing is present. If the force or torque acting on a stationary motor exceeds the 
static friction, the motor effect should be present. In their experiment Eshuis et al.[T2] 
observed that there is a net rotation of the asymmetric vanes that is only non-zero if 
the bath particles are sufficiently agitated. This is consistent with our analysis that 



gives an expression for the torque acting on a stationary rotor, equation (24). If this 
is larger than the static friction, the rotor starts to turn. If the rotor is immobile the 
granular temperature of the bath particles can be increased by increasing the vibration 
amplitude or frequency until the threshold is reached. 

Eshuis et al.[T2] also observed a convective motion that accompanies the motor 
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Figure 12. Dimensionless density of the bath particles (second layer) versus 
the dimensionless distance to the spherocylinder center r/Rs, with a bath density 
corresponding to a two-dimensional coverage of 0.18. The bath particle radius is 
R = 1.5Rs- The black and red curves correspond to homogeneous {a = 0.2) and 
heterogeneous (a+ — 0.8, = 0.2) sphere cylinders, respectively. The total length of 
the spherocylinder is L+Rg = 16i?s, The dashed line represents the Boltzmann-Lorentz 
assumption of a perfect homogeneous bath. 

rotation. This is different from the effect observed in a system of vertically vibrated 
glass beads where a toroidal convection roll with a vertical axis appears [27] . The origin of 
this effect is the inelastic collisions between particles and side wall [28]. In the experiment 
of Eshuis et al., however, convection accompanies the rotation of the vanes and is not 
the origin of the motor effect. 

A more quantitative analysis of the influence of dry friction on the ratchet effect 
requires its incorporation in the kinetic description. We are currently pursuing this 
direction. 

We thank J. Piasecki for useful discussions. 
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